A hollow sphere as a detector of gravitational radiation 
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The most important features of the proposed spherical gravitational wave detectors are closely 
linked with their symmetry. Hollow spheres share this property with solid ones, considered in the 
literature so far, and constitute an interesting alternative for the realization of an omnidirectional 
gravitational wave detector. In this paper we address the problem of how a hollow elastic sphere 
interacts with an incoming gravitational wave and find an analytical solution for its normal mode 
spectrum and response, as well as for its energy absorption cross sections. It appears that this 
shape can be designed having relatively low resonance frequencies (~ 200 Hz) yet keeping a large 
^ , cross section, so its frequency range overlaps with the projected large interferometers. We also 

' apply the obtained results to discuss the performance of a hollow sphere as a detector for a variety 

' of gravitational wave signals. 
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■ I. INTRODUCTION 

I I Thirty- five years after the beginning of the experimental search for cosmic gravitational waves (GW), several 
bJO- resonant-mass detectors (cryogenic cylindrical bars) are currently monitoring the strongest potential sources in our 
^ ' Galaxy and in the local group [Q. The sensitivity of such detectors is ft. ~ 6 x 10""'^^ for millisecond GW bursts, or, 
•rH , in spectral units, 10^^^ Hz^^/^ over a bandwidth of a few Hz around 1 kHz. A further improvement in sensitivity 
/\ ' and bandwidth is expected from the operation at ultralow temperatures of the two bar detectors NAUTILUS |^ 
^ , and AURIGA |^ in Italy, and even better sensitivities and bandwidths will come about as more advanced readout 
■ - - ' systems are developed. Projects for spherical resonant-mass GW detectors have emerged in the last few years in the 
resonant-mass community , due to their remarkable advantages with respect to the operating bars Q . 

In a cylindrical bar only the first longitudinal mode of vibration interacts strongly with the GW, and consequently 
only one wave parameter can be measured: the amplitude of a combination of the two polarization states. On the 
other hand each quadrupole mode of a spherical mass is five-fold degenerate (its angular dependence is described in 
terms of the five spherical harmonics Y/m(^,</') with 1=2 and m = — 2,...,2), and presents an isotropic cross section. 
The cross section of the lowest order (n=l) mode is the highest, and is larger than that of a cylindrical antenna made 
of the same material and with the same resonant frequency by a factor of about 0.8 {Rs/Rbf' 1^10 1 where Rs and Rb 
are the radius of the sphere and of the bar, respectively. This means a factor of 20 over present bars. Moreover, the 
sphere's cross section is also high at its second quadrupole harmonic. 

The five-fold degeneracy of the quadrupole modes enables the determination of the GW amplitudes of two polariza- 
tion states and the two angles of the source direction. The method first outlined by Forward |^ and later developed by 
Wagoner and Paik [|lO| , consists in measuring the sphere vibrations in at least five independent locations on the sphere 
surface so as to determine the vibration amplitude of each of the five degenerate modes. The Fourier components of 
the GW amplitudes at any quadrupole frequencies and the two angles defining the source direction can be obtained 
as suitable combinations of these five outputs l5P,|^, |ll| , |l2| . 

The signal deconvolution is based on the assumption that in the wave frame (that in which the z axis is aligned 
with the wave propagation direction) only the 1=2 and m=±2 modes are excited by the GW, as the helicity of a GW 
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is 2 in General Relativity. One can take advantage of this to deconvolve the wave propagation direction and the GW 
amplitudes in the wave frame. 

Most of the nice properties of a spherical GW detector depend on its being spherically symmetric. A spherical shell, 
or hollow sphere, obviously maintains that symmetry, so it can be considered an interesting alternative to the usual 
solid sphere. In order to have a good cross section, a resonant GW detector must be made of a high speed of sound 
material, and have a large mass. The actual construction of a massive spherical body may be technically difficult. In 
fact, fabricating a large hollow sphere is a different task than fabricating a solid one. Casting a hollow half sphere is a 
nearly two dimensional cast, at odds with casting a solid sphere, which requires rather special moulds. As an example 
of the feasibility of large two dimensional casting we can mention the fabrication of propellers of more than 10 meters 
in size and masses of the order of 100 tons |Q. Two hollow hemispheres could then be welded together with electron 
beam techniques. However, while it is known that these welding technique preserve most of the properties of the bare 
material, its effect on the acoustic quality factor (a relevant paramenter in resonant mass detectors) must be further 
studied. 

We have investigated the properties of a hollow sphere as a potential GW antenna. The purpose of this paper is 
to present a detailed report of the main results of such an investigation, and to discuss the real interest of this new 
detector shape. 

In section 2 we present the complete analytical solution of the eigenmode problem for a hollow sphere of arbitrary 
thickness, including the full frequency and amplitude spectrum. Section 3 is devoted to the cross section analysis, 
while in section 4 we take up the study of the system sensitivity to various GW signal classes. Finally, we present an 
outlook and summary of conclusions in section 5. 



II. NORMAL MODES OF VIBRATION AND EIGENFREQUENCIES OF A HOLLOW SPHERE 

In this section we consider the problem of a hollow elastic sphere in order to obtain its normal modes and frequency 
spectrum. This is a classical problem in Elasticity theory which was posed and partly addressed already in the last 
century, see e.g. Q and references therein. 

Let R and a be the outer and inner radius of the sphere, respectively. The elastic properties of the sphere, provided 
it is homogenous and isotropic, will be described by its Lame coefficients, A and /i, and its density, p. As is well known 
(see, e.g., [pi), the normal modes are obtained as the solutions to the eigenvalue equation 

V^u + (1 + A//x) V(V ■ u) = -fc^u, {k"^ =uj^p/p,) , (2.1) 

subject to the boundary conditions that the solid's surface be free of any tensions and/or tractions; these are 
expressed by the equations 

a,jnj=0 at r^R and at r = A (i? > a > 0), (2.2) 

where the sphere's surface S has outward normal n. The possibility of a spherical shell (a = i?), and that of a solid 
sphere (a — 0), are allowed. The stress tensor cry is given by Q 

aij ^\uk,kSij + 2pLU(i^j). (2.3) 



The general solution to (2.1) can be cast in the form 

u(x) = Co V(/)(x) + iCi LV^x) + iC2 V X LV^x) + 

Do V0(x) + iDi L?/^(x) + iZJa V X LV;(x), (2.4) 

where Ci, Di are constants, L = — ix x V, and the scalar functions cj), ip, (f), ^ are given by 

^{^) = jiiqr)Yi„,{9,ip), i^i^) = jiikr)Yi„,{9,ip), (2.5) 
0(x) = yi{qr)Yi^{9,^), ^^(x) - yi{kr)Yi^{9,^), (2.6) 



where q = k^J p/ {\ + /i), and Yim denotes a spherical harmonic. Finally, ji and yi are the standard Bessel functions 
of the first and second kind, respectively (see, e.g., [^). The latter (which are singul ar a t the origin) must be 
included in our case, as r = lies outside the boundary S. The boundary conditions (2^) become, after rather 
lengthy calculations, a system of linear equations which splits up into a 4 x 4 linear system for (Co, C2, Do, D2), and 
a 2 X 2 system for (Ci, Di). That is, we have a linear system of the form: 
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Ap 
At 



Cp 



= 0, 



(2.7) 



with 



Cp = (Co, C2,Do, D2Y, 



(Ci,Di 



where the superscript t denotes transpostion, and the corresponding matrices are: 



/ PiiqR) -l{l + l)s-^l3i{kK) PjiqR) -l{l + l)s-^f3i{kR)\ 



Ap = 



-s-^PsikR) (3i{qR) 
-l{l + l)s-^Piika) Mqa) 
-s-'^Psika) hiqa) 



-s^MkR) 
-l{l + l)s-^Pi{ka) 



(2.9) 



and 



At = 



[3i{kR) /3i{kR) 
(3i{ka) Piika) 



Here s = q/k, and we have introduced the set of functions: 



P2{Z)^J['{Z) 



Psiz) = i/32(z) - {1 - ^} /3o(^), f34{z) = p2{z) ~ ^z^Poiz), 



(2.10) 



(2.11) 
(2.12) 



while the tilded ones are their singular counterparts, with yi instead of ji (i.e., (3o{z) = yi{z)z~''^, and so on). The 
matrices Ap and At are functions of kR, and depend on the parameter a/R, and, in the case of Ap, also on s ^. 
The discrete set of kR values that make compatible the system (2.7) constitute the spectrum of the elastic sphere. 
We can distinguish two families of normal modes: 
(i) Toroidal modes. These are characterized by 

detAT = 0, Cp=0. (2.13) 

Hence they are purely tangential, and their frequencies depend only on the ratio a/R. Their amplitudes are 

uli^{^)^Tni{r)tLYUO,cb), (2.14) 



with 



Tniir) = Ci{n,l) |/3i(fc„;i?)j;(A;„;r) - /3i(fc„;i?)y;(fc„;r)| 



(2.15) 



where C'i{n,l) is fixed b y the chosen normalization. The corresponding eigenvalues are obtained as solutions to the 
trascendental equation ( 2.13 ). For the degenerate limit a = R the equation to be solved is 



, /3i(fci?) P,{kR) \ 
(3[{kR) P[{kR) ) 



(2.16) 



with the prime denoting differentation respect to the argument. Using standard properties of Besscl functions p5[ , it 
can be easily shown that 

(3i{kR)i3[{kR) - Pi{kR)(3[{kR) = {kR)-%kRf + 2-1(1 + 1)], 

and, in this case, there is only one eigenvalue for each I > 1, given by the only root of the above equation, (kiR)^ = 
1(1 + 1) — 2^ Figure |l| displays fc„;i? as a function of a/R for the first few toroidal modes. The existence of just 
one mode for each Z > 1 in the thin shell limit shows as a divergence of kniR when a/R approaches 1 and n > 1. In 
figure 1^ we plot the normalized toroidal amplitudes Tni (r) for two quadrupolar modes and three different values of 
the parameter a/R. We observe that their absolute values at the outer surface show little dependence on the ratio 
a/R. 
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FIG. 1. Functional dependence of the first few toroidal eigenvalues of a hollow sphere on the ratio a/R. Solid sphere values 
(a = 0) are found on the intersections with the ordinate axis. 
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FIG. 2. Toroidal mode radial functions for the first two quadrupole harmonics and a few values of the geometric ratio a/R. 
The magnitude represented in abscissas is such that the region plotted spans radially the material thickness of the hollow 
sphere. 
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(ii) Spheroidal modes. This second family is charaterized by: 



dct Ap = 0, Ct = 0. 



(2.17) 



In this case, the expressions get mor e inv olved, as we have to handle a 4 x 4 determinant. Once the spectrum fc„; is 
found for given a/R and s, the system ( [2.7| ) can be solved for C2/C0, Dq/Cq, and £'2/^0. If we label these coefficients 
Po{n,l), pi{n,l), p2{n,l), the eigenmodes can be written as 



= Nni{r)Yhr,{e, (P)n - iE^i{r)n x Ly,„, (61, . 



(2.18) 



with 



iV„,(r) = C„(n,() 



qnir 



Pi y'lilnir) - p2 [n, l)^-^—^^ Vi (knir) 



Eni{r) = Coin, I) [jiiqnir) ~ Po{n,l){knir ji{knir)}' - 

Pi{n,l)yi{q„ir) - P2{n,l){knir yi{knir)}'] , 



(2.19) 



(2.20) 



where Co{n, I) is, again, free up to normalization. The spectrum for the degenerate case a = i? is given by the solutions 
to 



( (3i{qR) -l{l + l)s-^Pi{kR) PjiqR) -1(1 + l)s-^i3i{kR) \ 



det 



piiqR) 



^PsikR) PiiqR) 



-s^MkR) 



PiiqR) -lil + l)s-^P[ikR) P'iiqR) -lil + l)s-^f3[ikR) 



XP'iiqR) 



-0, 



(2.21) 



-'P'sikR) I 



which happens to have two solutions for each value of I when / > 1 and only one root for / < 2 ^ 

Plotting kniR as a function of a/i?, we see that the third and higher roots diverge as the inner radius approaches i?, 
see figures || and ^ Figures show the normalized radial functions for a few spheroidal modes and values of ajR. 
As in the toroidal case, their values at r = i? (where measurements using transducers are to be made eventually) are 
nearly independent of ajR. 



III. CROSS SECTION FOR THE HOLLOW SPHERE 



A convenient way to characterise a resonant detector sensitivity is through its GW energy absorption cross section, 
defined as 

where Ai?Q(a;) is the energy absorbed by the detector at frequency and ^iio) is the incident flux density expressed 
e.g. in watt/m^ Hz. Estimation of crabs('^) requires a hypothesis about the underlying gravitation theory to calculate 
$(0;), and specification of the antenna's geometry to calculate AiJ(j((jj). Here we shall assume that General Relativity 



^This parameter is a function of the Poisson ratio a; for the usual value a = 1/3, s — 0.5 and \/ ^ = 2. These values are 
assumed, unless otherwise stated. 

^ This equation shows explicitly a property shared by all toroidal modes, namely, that their dimensionless eigenvalues k^iR 
do not depend on the elastic properties of the material. 

^The purely radial case Z = is simpler, because the eigenvalue equation (protect 2.21) becomes 

P4{qR)Pi{qR) ~ A{qR)i3iiqR) = 0, 



and has only one solution, namely qR = (/i/A)-^3 — ^/A. Unlike toroidal eigenvalues, spheroidal ones do depend on /i/A. 
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FIG. 5. Spheroidal mode radial A^-functions — see equation (2.19) — for the first two monopole harmonics and a few values 
of the geometric ratio a/ R. The magnitude represented in abscissas is such that the region plotted spans radially the material 
thickness of the hollow sphere. 
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FIG. 6. Spheroidal mode radial functions — see equations (2.19) and (p.20|) — for the first quadrupole harmonic 
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FIG. 7. Spheroidal mode radial functions — see equations (2.19) and (2.20) — for the second quadrupole harmonic. 



is the correct gravitation theory, and proceed to calculate the oscillation energy of the solid as a consequence of its 
excitation by an incoming GW, which we shall naturally identify with AEa{uf). We briefly sketch the details of the 
process now. 

As shown in an elastic solid's response to a GW force can be expressed by a very general formula, which is 
easily particularised to a spherically symmetric body such as the solid sphere or the hollow sphere. In both cases, as 
we have just seen, the vibration eigenmodes belong into two families (spheroidal and toroidal), but GWs only couple 
to quadrupole spheroidal harmonics. If the frequencies of these modes are noted by a;„2 {n= \ for the lowest value, 
n = 1 for the next, etc.) and the corresponding wavefunctions by Vin2m 

(x) then the elastic displacements are given by 



u(x,t) = ^ 



U"2m(x)gi"(0 



m=-2 



where 



(3.2) 



g''"'Ht') smujn2{t-t')dt' 



(m = -2,...,2) 



(3.3) 



and g'''^\t) are the quadrupole components of the Riemann tensor, while &„ is an overlapping integral factor of the 
GW's tidal coefficient over the solid's extension. Much like in the case of a solid sphere, it has dimensions of length, 
and is given by a definite integral of the radial terms in the wavefunction u„2m (x) ; more specifically. 



R 



P_ 

M 



" [Nn2 (r) + 3i?„2 (r)] dr = - ^^^^ [G2 {R) - G2 (a)] , 



47r q„2i? 



where we have introduced the dimensionless function 
.3 



G2(z) = 



[j2{qn2z) + pi{n,2)y2{qn2z) 3po{n,2)j2{kn2z) - 3p2in,2)y2{kn2z)] 



i?3 _a3 

and assumed the following normalization for the wavefunctions 



Solid 



|u„,,„|' pd'^= / r'drp [N?^ir)+l{l + l)Elir)] =M 



(3.4) 



(3.5) 



(3.6) 



The calculation of Ai?a(^) can now be pursued along the lines set up in reference the Fourier transform U(x, uj) 
of the response function u(x, t) of equation ( |3.2| ) is calculated, whereby the spectral energy density can be obtained as 
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FIG. 8. Cross sections of a hollow sphere in its first two quadrupole modes as a function of thickness. Values are referred to 
the cross section of a solid sphere in its first quadrupole resonance, whose radius is assumed to be equal to the outer radius of 
the hollow sphere. 



W{uj) 



1 

T 



Solid 



icj2 |U(x,w)|^ pcfx 



(3.7) 



where T is the integration time of the signal in the detector. The energy deposited by the GW in the n-th quadrupole 
mode is hence calculated by integration of this spectral density over the linewidth of the mode. It is readily found 
that 



1 ^ I 



(3.8) 



where G(™)(^) is 

the Fourier transfor m of g'^'^'^jt). 

The GW flux in the denominator of ( |3.l| ) is (clearly) proportional to the sum in the rhs of (3.S), the proportionality 
factor being in turn proportional to uf^ — see for a detailed discussion — , so we finally obtain 



0'abs('^n 



15 



(3.9) 



where = ^/p, M is the detector's mass, and G is Newton's constant. This equation allows relatively easy numerical 
evaluation of the cross sections, as well defined computer programmes can be written for the purpose. 

As we have seen in section 2 above, the eigenval ues a nd wavefunctions of a hollow sphere only depend on the ratio 
a/R, and therefore so does the quantity (fc„2&n) in (3^). So the cross section cr„ only depends on that ratio, too, once 
a suitable unit of mass is adopted for reference. In figures^ and ^ we plot (t„ for the first two quadrupole modes of 
the hollow sphere in two different circumstances: in figure Wwe assume a hollow sphere of fixed outer radius — thus its 
mass decreases with thickness — , and in figure ^ we have instead assumed that the mass of the hollow sphere is fixed, 
so that its geometrical size increases as it gets thinner. In either case we see that, for the higher mode, the maximum 
cross section does not happen at a=0, but at some intermediate inner radius: for a « 0.37745i?, the cross-section for 
the second quadrupole mode equals that of the first, and we have the possibility of working with a detector with the 
same (high) sensitivity at two frequencies. 



IV. SENSITIVITY TO GW SIGNALS 



We assume that the mechanical oscillations induced in a resonant mass by the interaction with the GW are 
transformed into electrical signals by a set of identical noiseless transducers (for the sake of simplicity, we consider 
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here non-resonant transducers), perfectly matched to electronic amplifiers with noise temperature T„. Unavoidably, 
Brownian motion noise associated with dissipation in the antenna and electronic noise from the amplifiers limit the 
sensitivity of the detector. We refer the reader to Jl6|-p^ for a complete discussion on the sensitivity of resonant-mass 
detectors and report here only a few basic formulas for the evaluation of the detector sensitivity to various signals. 

The total noise at the output of each resonant mode can be seen as due to an input noise generator having spectral 
density of strain Sh{f), acting on a noiseless oscillator. Sh{f) represents the input GW spectrum that would produce 
a signal equal to the noise spectrum actually observed at the output of the detector instrumentation. In a resonant- 
mass detector, this function is a resonant curve and can be characterized by its value at resonance Sh{fn) and by its 
half height width. Sh{fn) can be written as: 

Here Te is the thermodynamic temperature of the detector plus a back-action contribution from the amplifiers, and 
Qn is the quality factor of the mode. 

The half height width of Sh{f) gives the bandwidth of the resonant mode: 

A/„ = ^r-V2 (4.2) 
Here, r„ is the ratio of the wideband noise in the n-th resonance bandwidth to the narrowband noise, 

Tn 
'^f^n Qn 

where /9„ is the transducer coupling factor, defined as the fraction of the total mode energy available at the 
transducer output. 

In practice r„ <^ 1 and the bandwidth is much larger than the pure resonance linewidth fn/Qn- In the limit r„ 0, 
the bandwidth becomes infinite. The bandwidth of the present resonant bars is of the order of a few Hz If a 
quantum limited readout s^y^tcm were available, values of the order of 100 Hz could be reached |19|j20||. 



The equations ( 4.1 and ]L2) can be used to characterize the sensitivity of the quadrupole modes of a hollow spherical 
resonant-mass detector. The optimum performance is obtained by filtering the output with a filter matched to the 
signal. The energy signal-to-noise ratio (SNR) of the filter output is given by the well-known formula 

/ + 00 I TTf r\l2 
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TABLE I. Main features and sensitivities for several hypothetical hollow spheres of two different materials. 



M (ton) 2R{m) (7? - a) (cm) /i (Hz) /2 (Hz) cti (m^ Hz) aa (m^ Hz) ^/S^(Hz-^/^ 



CuAl 



A15056 



200 
200 
100 
100 
40 
40 
200 
100 
40 
40 



81 
25 
31 
12 
22 
11 
90 
37 
35 
14 



395 
191 
302 
185 
399 
281 
273 
230 
361 
218 



1188 
753 
1161 
738 
1543 
1115 
935 
896 
1370 
866 



1.5 10"^^ 

1.1 10"^^ 
5.8 10-^* 

5.6 10"^* 
2.3 10"^* 

2.2 10"^* 
1.8 10"^^ 

7.7 10"^* 
3.2 10"^* 
3.0 10"^* 



2.4 10"^'' 

1.7 10^2^ 

8.8 10"^" 

8.2 10"^" 

3.5 10"^'' 

3.3 10"^'' 

2.9 10~^^ 
1.2 10"^^ 
4.8 10"^" 

4.4 10"^" 



4.7 10"^" 

5.4 10^^'' 

7.5 10"^'' 
7.7 10"^" 
1.2 10"^^ 

1.2 10"^^ 

4.3 10"^" 

6.6 10"^'* 
1.0 10"^^ 
1.0 10"^^ 



3.7 10"^* 
4.4 10~^* 
6.1 10"^* 

6.3 10"^* 
9.7 10"^" 
1.0 10"^^ 

3.4 10"^* 
5.4 10"^* 
8.3 10"^* 
8.7 10"^* 



where H{f) is the Fourier transform of h{t). 

We now report the SNR of a hollow spherical detector for various GW signals. To be specific, we shall assume that 
the thermodynamic temperature of the detector can be reduced to below 50 mK, and that the quality factors of the 
modes are of the order of 10^, so that the overall detector noise will be dominated by the electronic amplifier noise. 
If we express the energy of the latter as a multiple of the quantum limit, i.e., kTn = NTiuj then the strain spectral 
density becomes 



G 4nf3nNh 

^h\Jn) — 



(4.5) 



In these conditions the fractional bandwidth /S.fn/ fn becomes of the order of /3„ that we assume of about 0.1. We 
shall consider hollow spheres made of the usual aluminium alloy Al 5056 and of a recently investigated copper alloy 
(CuAl) . Table | displays numerical values of the most relevant parameters for a few example detectors with a 
noise level equal to the quantum limit, i.e. = 1. 



A. Bursts 



We model the burst signal as a featureless waveform, rising quickly to an amplitude /iq and lasting for a time Tg 
much shorter than the detector integration time Ai = A/^^. Its Fourier transform will be considered constant within 
the detector bandwidth: H{f) ~ H{fn) = Hq. From (O) we get 



SNR = 



2TTAfnH^ 

Shifn) 



For SNR = 1, and using the equation Hq^ 



/iS"'Vg, we find 



Shifn) 



-,1/2 



2^A/„ 



(4.6) 



(4.7) 



The level /ij^™ 



10~^^ can be reached by the lowest order mode of a typical large hollow spherical detector 
such as the one being considered. The GW luminosity of burst sources is still largely unknown, so it is difficult to 
accurately estimate their detectability. The above sensitivity is however likely to enable the detection of GW collapses 
in the Virgo cluster for an energy conversion of lO"** M© into a millisecond GW burst. See table || for a few specific 
examples. 



B. Monochromatic signals 



We consider a sinusoidal wave of amplitude ho and frequency fs constant over the observation time t„ 
transform amplitude at /„ is ■^hotm with a bandwidth given by t^. The SNR can be written as 



The Fourier 
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TABLE II. Sensitivity to burst and monochromatic (integrated for one year) GW signals of a few hollow spheres of two 



different materials. 





M (ton) 


27? (m) 


{R — a) (cm) 


,fi (Hz) 


.f2 (Hz) 


\ "1 /burst 


\ "2 /burst 




ru 




m 


CuAl 


200 


6 


25 


191 


753 


4.9 10"^^ 


2.0 10-2^ 


1.4 10" 


-27 


1.1 10" 


-27 




100 


4 


31 


302 


1161 


5.5 10-22 


2.3 10-22 


1.9 10" 


-27 


1.5 10" 


-27 


A15056 


100 


6 


37 


230 


896 


5.5 10-22 


2.3 10-22 


1.7 10" 


-27 


1.4 10" 


-27 



SNR 



{i + r„ [Q^ (1 - + /.V/^] } 



For SNR — 1 we obtain a minimum detectable value of ho, which at fs = fn is 



1/2 



(4.8) 



(4.9) 



See table || for a few specific examples. For instance, the nearby pulsar ||2^ PSR J0437-4715, at a distance of 150 
pc, might emit at 347 Hz a GW amplitude (optimistically) of 2 x IQ-^^. This would give SNR = 100 on a hollow 
spherical detector having M = 100 tons after integrating the signal for 1 year. 



C. Chirps 



We consider here the interaction of the hollow spherical detector with the waveform emitted by a binary system, 
consisting of either neutron stars or black holes, in the inspiral phase. The system, in the Newtonian regime, has a 
clean analytic behaviour, and emits a waveform of increasing amplitude and frequency that can sweep up to the kHz 
range of frequency. 

From the resonant-mass detector viewpoint, the chirp signal can be treated as a transient GW, depositing energy 
in a time-scale short with respect to the detector damping time |2^. We can then use (4.6) to evaluate the SNR, 
where the Fourier transform H{ fn) at the resonant frequency /„ can be explicitly written as 



H{fn) 



n 2 



h{t) cos(27r/„t) dt 



n 2 



h{t) sin(27r/„t) dt 



1/2 



(4.10) 



with h{t) indicating h^{t) or h^{t). Substituting into ( 4.10 ) the well-known chirp waveforms for an optimally 
oriented orbit of zero eccentricity in the Newtonian approximation pq| , the SNR for chirp detection is [p4[ : 



SNR 



21/35 ,2^/„ 1 



(4.11) 



Mc is the chirp mass defined as Mc — (mim2) ' (m-i + TO2)- ' , where mi and 7712 are the masses of the two 
compact objects and r is the distance to the source. The chirp mass is the only parameter that determines the 
frequency sweep rate of the chirp signal in the Newtonian approximation, and can be determined by a double passage 
technique |^^: much like in a solid sphere detector, one can measure the time delay T2 — ti between excitations of 
the first and second quadrupole modes on a hollow spherical detector to calculate the chirp mass through equation 



M, = 28/5 ( — - 
V256 / G 



-8/3 



,-8/3 ~ 



3/5 



T2 - Tl 



(4.12) 



where wi and uj2 are the angular frequencies of the first and second quadrupole modes, respectively. Time delays 
are of the order of a fraction of a second for the hollow spheres considered in this paper, well within the timing 
possibilities of resonant mass detectors 

Another consequence of the multimode and multifrequency nature of a spherically symmetric detector is the pos- 
sibility to determine the orbit orientation by the measurement of the relative proportion of the two polarisation 
amplitudes, and thereby the distance to the source and the intrinsic GW amplitudes ||2^. See figs. |o| and [ill for a 
specific example referring to optimally oriented circular orbits. 
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FIG. 10. Contours of constant chirp mass Mc in mi, m2 space. At each chirp mass corresponds the maximum distance r at 
which the chirp can be observed with a SNR=10 by a 200 ton CuAl hollow sphere, 6 meters in diameter, at its first resonance 
frequency /i = 191 Hz. The reported chirp mass values (in units of solar masses) and the corresponding maximum distances 
are: a : = 8.0, r = 214 Mpc, h : Mc = 4.0, r = 119 Mpc, c : Mc = 2.6, r = 84 Mpc, d : Mc = 1.2, r = 45 Mpc. 



Because of the Newtonian approximation, eqs. (4.11) and (4.12) become inaccurate near coalescence. In analogy 
with previous analyses ^3 24|, we limit our considerations to the frequency at which there are still five cycles remaining 
in the waveform until coalescence. The highest chirp mass values reported in the figures are determined by the 
requirement that the five-cycle-frequency of the source be larger than the resonant frequencies of the detector. 



D. Stochcistic background 

In this case h{t) is a random function and we assume that its power spectrum, indicated by Sgwif), is fiat and its 
energy density per unit logarithmic frequency is a fraction flgw{f) of the closure density pc of the Universe: 



dpgw 



dlnf 

Sgwif) is given by 

2G 



= ^gwPc (4-13) 



SgM) = —r'^sM)Pc (4.14) 



TT 



The measured noise spectrum Sh{f) of a single resonant-mass detector automatically gives an upper limit to Sgwif) 
(and hence to flgwif))- 

Two different detectors with overlapping bandwidth A/ will respond to the background in a correlated way. The 
SNR of a GW background in a cross correlation experiment between two detectors located near one another and 
having a power spectral density of noise Sj^{f) and Sf^{f) is p^ ]: 



2 \ 1/4 



SNR = I #^A/t,„ ) (4.15) 



where tm is the total measuring time. 

Detectors located some distance apart do not correlate quite so well because GWs coming from within a certain 
cone about the line joining the detectors will reach one of them before the other. The fall-off in the correlation with 
separation is a function of the ratio of the wavelength to the separation, and has been studied for pairs of bars, pairs 
of interferometers P^'la ^^'^ P^i^ spherical detectors [M . 



13 



10 



a 




FIG. 11. Contours of constant Mc for the same hollow sphere as in fig. y_0[ observing the chirp at the second resonance 
frequency /2 — 753 Hz, with SNR=10. The chirp masses and the maximum distances are: a : Mc = 2.0, r = 131 Mpc, 
b : Mc — 1.2, r = 86 Mpc, c : Mc = 0.9, r — 65 Mpc. If the double passage technique is applied, the delay times between the 
excitation of the first and the second mode by the chirps of the given mass are: a : 160 ms, b : 373 ms, c : 648 ms. 



Assuming two identical large hollow spherical detectors are co-located for optimum correlation, the background will 

1/2 /in7 \ 1/2 

^ 10 sec ^ 



reach a SNR = 1 if flgw is 



-9 f fn \^ ( V Slifn) 



20 Hz 



[200 Hz j I 10-24 HZ-V2 I I 10-24 Hz-1/2 / V A/„ 



(4.16) 



where the Hubble constant has been assumed 100 kms~^. 

Hollow spherical detectors can set very interesting limits on the GW background. In particular, following recent 
estimations based on cosmological string models [Q, it emerges that experimental measurements performed at the 
level of s ensit ivity attainable with these detectors would be true tests of Planck-scale physics. 

Eqs. (4.15) and ( 1.16| ) hold for whichever cross-correlation experiment between two GW detectors adjacent and 
aligned for optimum correlation. An interesting consequence is that the sensitivity of a hollow sphere-interferometer 
observatory will be unprecedented. It can be worth to build a hollow spherical mass detector close to a large 
interferometer, like LIGO or VIRGO, to perform stochastic searches [pll]. 



V. CONCLUSIONS 



In this paper we have been mainly concerned with the problem of how an elastic hollow sphere responds to a GW 
signal impinging on it. To address this problem we have developed an analytical procedure to fully sort out the 
eigenfrequencies and eigenmodes of that kind of solid, then applied it to calculate the GW absorption cross section 
for arbitrary thicknesses and materials of our solid. 

When realistic hypotheses are made regarding the size and material of a possible GW detector of this shape, we 
have seen that a hollow sphere can be advantageous in several respects. It has all the features associated with its 
symmetry, such as omnidirectionality and capability to determine the source direction and wave polarisation. Also, 
its quadrupole frequencies are below those of an equally massive solid sphere, thus making the low frequency range 
accessible to this antenna with good sensitivity. We have investigated the system response to the classical GW signal 
sources (bursts, chirps, continuous and stochastic) for several sizes and materials, and seen that interesting signal-to- 
noise ratios are attainable with such a detector. Also, its bandwidth partly overlaps with that of the projected large 
interferometers ||3^,Q , so potentially both kinds of detectors can be operated simultaneously to make hybrid GW 
observatories of unprecedented sensitivity and signal characterisation power. 
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While it seems possible to cool a 100 ton solid sphere down to 50 mK ||3J], the possibility of cooling a large 
hollow sphere at such low temperatures, as well as the fabrication technique and the influence of cosmic rays on a 
low-temperature GW detector of that shape and dimensions, are currently under investigation. 
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